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Abstract. We explore the relationship between the nonlinear matter power spectrum and 
the various Lagrangian and Standard Perturbation Theories (LPT and SPT). We first look at 
it in the context of one dimensional (1-d) dynamics, where ILPT is exact at the perturbative 
level and one can exactly resum the SPT series into the ILPT power spectrum. Shell crossings 
lead to non-perturbative effects, and the PT ignorance can be quantified in terms of their 
ratio, which is also the transfer function squared in the absence of stochasticity. At the 
order of PT we work, this parametrization is equivalent to the results of effective field theory 
(EFT), and can thus be expanded in terms of the same parameters. We find that its radius of 
convergence is larger than the SPT loop expansion. The same EFT parametrization applies 
to all SPT loop terms and if stochasticity can be ignored, to all N-point correlators. In 3- 
d, the LPT structure is considerably more complicated, and we find that LPT models with 
parametrization motivated by the EFT exhibit running with k and that SPT is generally a 
better choice. Since these transfer function expansions contain free parameters that change 
with cosmological model their usefulness for broadband power is unclear. For this reason 
we test the predictions of these models on baryonic acoustic oscillations (BAO) and other 
primordial oscillations, including string monodromy models, for which we ran a series of 
simulations with and without oscillations. Most models are successful in predicting oscillations 
beyond their corresponding PT versions, confirming the basic validity of the model. We show 
that if primordial oscillations are localized to a scale q, the wiggles in power spectrum are 
approximately suppressed as exp[—A;^S^(g)/2], where Ti{q) is rms displacement of particles 
separated by q, which saturates on large scales, and decreases as q is reduced. No oscillatory 
features survive past k ~ 0.5h/Mpc at z = 0. 
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1 Introduction 

Effective Field Theory (EFT) approach to large scale structure (LSS) [1-7] has lately received 
a lot of attention as a way to extend the validity of cosmological perturbation theory (PT) 
(see e.g. [8-10]), at a cost of introducing free parameters (e.g. [2, 5, 6]). These parameters 
must obey some symmetry requirements. For example, any local nonlinear scrambling of 
matter must obey mass and momentum conservation [11], and one can show that at lowest k 
(where k is the wavevector amplitude of the Fourier modes), the leading order effects scale as 
ak'^Plik), where Phik) is the linear matter power spectrum. In the halo model [12-15] one 
can assemble the local scrambling of dark matter into dark matter halos, and one can perform 
a Taylor series expansion of the halo profile to extend this into a series of even powers of k 
[16]. There is a second term that in EFT language is called the stochastic or mode coupling 
term, for which mass and momentum conservation require to initially scale as k^. This term 
is often ignored in EFT calculations, but it should be the dominant term on small scales. 
Since PT calculations already enforce the conservation laws such terms can be applied to 
any PT scheme as a way to correct for whatever is missing in PT calculations. So EFT can 
be viewed simply as a parametrization of the ignorance of a given PT model, absorbing any 
discrepancy between the true solution and the PT solution, and parametrizing it in terms of 
a simple parameter expansion. 

The usefulness of this approach then depends on the convergence radius of this expansion, 
or more simply, the range of k over which the lowest order EFT term(s) restore the exact 
solution. The answer will depend on the specific PT model implementation: one can apply 
PT ignorance to linear theory, for example, but that will not be very useful and will lead to 
a strong running of the EFT parameter with scale, as well as a large stochasticity [17, 18]. 
This is because the lowest order EFT correction is of the same order as PT at the 1-loop 
order, so it makes sense to add lowest order EFT term to 1-loop PT. Moreover, we will 
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argue that the convergence radius of EFT expansion should be a bit larger than that of loop 
expansion. While this discussion is general, a very relevant question is which PT to use. In 
literature we have PT approaches both in Eulerian space (SPT) (see e.g. [6, 8, 10, 19], and 
in Lagrangian space (LPT) (see e.g. [5, 9]), and at several different orders. These give very 
different predictions for the power spectrum, and it is unclear which is more successful. The 
purpose of this paper is twofold. One is simply to test the various PT models, parametrize 
their ignorance against the true answer in terms of EFT parameters, and study their scale 
dependence. The less scale dependence there is, the more useful the expansion. Our goal is 
to test several PT models, including some introduced here for the hrst time. 

Moreover, scale dependence alone is not the only criterion, as it could be a coincidence 
that the EFT parameter is roughly constant over a certain range of k. Our second purpose is 
to test the validity of EFT+PT approaches by applying it to modeling of oscillatory features in 
linear power spectrum. Baryonic acoustic oscillations (BAO) are a prime example of usefulness 
of EFT approach: BAO appear at low k, where we would expect EFT corrections to be valid. 
Since EFT corrections at the lowest order scale with PL{k), and PL{k) contains BAO, one 
would expect EFT correction to carry the signature of BAO. EFT correction is a small effect 
on top of a small BAO effect, so we would not expect it to be visible in simulations at low k, 
where sampling variance errors dominate. For this reason we ran a series of simulations with 
and without oscillations, but with the same initial conditions (as described in [20]), so that 
the sampling variance errors cancel [21]. 

In this paper, we adapt the simple parametrization of the dark matter overdensity via 
the density transfer function (see e.g. [17]). In this approach, EFT parameters are obtained 
by expanding the transfer function in even powers of k. This equivalence of EFT and transfer 
function expansions holds at the lowest order in PT (see e.g. [6]), but not beyond that. This 
means that the leading coefficient of transfer function expansion is A:^, but beyond that one 
can have an arbitrary Taylor expansion. Nevertheless, for simplicity we will try even powers 
of k as our expansion basis. Alternatively, one could also try different transfer functions 
at higher orders in density (this was recently indicated in [18]). With all this in mind, for 
simplicity we will refer to all such expansions as the EFT expansions, and parameters we will 
call EFT parameters. 

The outline of the paper is as follows. The relationship between LPT, SPT, and EFT is 
particularly simple in 1-d dark matter dynamics and in section 2 we hrst analyze this example 
using results from [22]. We then move to 3-d analysis in section 3, extracting EFT parameters 
for various PT models. In section 4 we apply these to the modeling of BAO features and 
other oscillations, and in section 5 we discuss the general lessons for the modeling of power 
spectrum. We present the conclusions in section 6. In appendix A we present the details of 
the construction of no-wiggle power spectrum. 

2 Investigation of the 1-d example 

To introduce the different schemes and their relation we will hrst look at the 1-d example. It 
is useful to look at the PT expansions in the context of 1-d dynamics, where ILPT (Zeldovich) 
solution is exact at the perturbative level, and can be shown to be identical to SPT in the 
inhnite loop limit [22], 


-PiLPT(fe, z) = '^ D+^~^^\z)PsPT,i-loop{k, z = 0), 
i=0 
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where i^sPT,i-ioop is the i-loop SPT term with 0-th loop given by PL(k), the linear power 
spectrum. These terms are multiplied with the appropriate growth rate D+{z), which we 
will for convenience normalize to D^(z = 0) = 1, so they can be dropped as long as we 
work at z = 0. From results of [22] we estimate the radius of convergence of SPT to be 
around (0.2-0.3)h/Mpc for ACDM like linear power spectra (in terms of power per mode): 
the scale at which the 10 loop SPT is accurate is /c ~ 0.2h/Mpc, and SPT series becomes 
highly oscillatory for k > 0.2h/Mpc. 

ILPT solution is however not the true solution to the dark matter Pdra{k), a consequence 
of shell crossings, at which ILPT sheets continue to stream through the shell crossings unper¬ 
turbed, while in the actual dark matter dynamics they stick together inside the high density 
sheets, the equivalent of halos in 3-d. In 1-d the ILPT solution is exact up to shell crossings, 
and it agrees with N-body simulations away from high density regions: the only difference 
is in the regions of shell crossings, which are broader in ILPT solution and typically have 
double peaks: ILPT solution artihcially spreads out the true density held. In the halo model 
language [12-14] the leading term correction at low k comes from applying this smearing to 
the true density prohles of dark matter, which gives rise to the even powers of k series cor¬ 
rections multiplying the linear power spectrum Pi^k), and the same is found in the context 
of EFT (e.g. [2, 5]). The leading order correction is thus the 2-halo (i.e. EFT like) ak'^Pi, 
with a positive sign to compensate for excessive smearing of ILPT displacements. Here 
corresponds to a typical scale of the streaming beyond the shell crossings, which can be sev¬ 
eral Mpc/h at z = 0. These stream crossing induced nonlinear corrections to PuPTik) are 
very large, 10% at k = 0.07h/Mpc and growing to a factor of 2 at fc = 0.3h/Mpc at z = 0. 
In 1-d it is clear that PT cannot address these stream crossings. Note however that all the 
deviations from ILPT are within a few Mpc/h. Hence, while the nonlinear effects are large 
in the power spectrum down to very low k, the correlation function can still be very close to 
ILPT on scales larger than a few Mpc (away from BAO peak) [22]. 

Even in the presence of stream crossings the ILPT held is well correlated with the dark 
matter held, a consequence of the fact that it spreads the high density peaks into a double 
peaked structure with an approximately constant radius independent of the position (and 
independent of the collapsed mass in the sheets). This can be quantihed by introducing the 
transfer function (see e.g. [17]) 


Pilpt(^) 


(i^lLPT^d m) 

((^ilpt^ilpt) ’ 


( 2 . 2 ) 


where (5ilpt and hdm are the ILPT and dark matter density perturbations in Fourier space, 
respectively. We can also introduce the cross-correlation coefficient 

2 ^ (<^lLPT<^dm)^ 

(hiLPT<^lLPT)(<^dm<^dm) 


It is expected to approach unity at low k and it has been shown that up to A: = 0.2/Mpc, 
the stochasticity 1 — t^lpt below 1% in 1-d [22]. We will denote auto-power spectrum as 
(27r)^h'^(k-|- ]i')Px{k) = (hx(k)(5x(k')), where X stands for dm, ILPT etc. and 5^ is Dirac 
delta function. The corresponding stochastic power Pj{k) is dehned as 


Pj{k) — [l “ ^iLPT(^)]-f’dm(A:). 


(2.4) 


The transfer function has to obey some symmetry properties, and in particular has to start 
as k'^ [11]. We will expand it into a general function of even powers of k. If r'iLPT(A:) < 1 then 
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Figure 1. Error of various models of 1-d power spectrum shown relative to the nonlinear simulations 
results. Show are PT results in black: ILPT/Zeldovich (dashed line), 1-loop SPT (dot-dashed line), 
2-loop SPT (double dot-dashed line), 5-loop SPT (long-dashed line) and linear theory result (dotted 
line). In addition we apply the transfer functions to these results giving us EFT+SPT and EFT+ILPT 
models in 1-d. Results for three different transfer functions are shown: going up to ai (in blue), 

(in red) and (in orange) in expansion given by Eq. (2.5). Thin grey horizontal dotted and dashed 
lines represent respectively 1% and 2% errors. Thin grey vertical solid lines represent maximal k 
values up to where EFT+ILPT models acheve 1% errors. Results are shown at redshift z = 0. 


1-d power spectrum 



we need a separate function riLPT(^) to fully describe the dark matter power spectrum given 
-Pilpt(^)) and if the goal is to specify the non-perturbative effects on the power spectrum 
then it is simpler to define, 

p = ^ilpt(^) = 1 + aiLPT(fe)fc^ ^ “*■ X] • (2-5) 

-PlLPT(fc) V ^ J 

Note that T{k) includes stochasticity and there is no guarantee that it can be expanded in 
terms of even powers of k, although we expect that at low k the leading term is aipuPT^^- 
We have fitted this expansion to the numerical results for Pdm(^) and Pilpt(^) given 
in [22]. This gives the values of the first three coefficients aipppT = 14(Mpc/h)^, a 2 ,iLPT = 
—4(Mpc/h)^ and Us^lpt = —40(Mpc/h)®. This is shown in top of figure 1 and the first 
coefficient is a good fit at 1% level up to ~ 0.2h/Mpc, while with three coefficients the fit 
is good to A: ~ 0.4h/Mpc. Note that we do not see the value aipLpr to approach a constant 
at low k: there is no sampling variance scatter, but there could be numerical issues with the 
simulations that prevent us from extracting the true value at very low k. We see the same 
issue in 3-d (recently also shown in [18]), and similar results were found for the displacement 
analysis of [23]. At very low k the nonlinear effects are really small, and these issues are 
unlikely to be relevant for any observations, since sampling variance errors are large on large 
scales. 

From figure 1 we see that we need just one EFT parameter even when go to 5th order 
in SPT (and beyond). This suggests that the radius of convergence in k for SPT expansion is 
smaller than for EFT expansion. Physically this makes sense: SPT breaks down for ~ 1, 
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and after that the collapse and shell crossings occnr (the same happens in 3-d, where one 
nsnally defines 5c ~ 1-68 as the linear density at collapse). There is some jnstification for the 
connting of EFT orders to be the same as the SPT loop orders, bnt the two can be different, 
and one expects EFT expansion to need fewer terms for a given SPT order. Finally, note that 
the stochastic term becomes of order 1% aronnd k ~ 0.2/Mpc at z = 0 and rapidly increases 
for higher k. In this regime the expansion in terms of even powers of k in eqnation 2 is not 
well jnstified, and instead one needs all powers of k, as in any Taylor expansion. 

To see the interplay between the SPT and EFT expansions, we can assnme oi^lpt 
applied to each SPT term and to fnll ILPT, or we order the terms assnming even powers of 
k of EFT expansion corresponds to a given SPT loop order. We can thns write, assnming 
?’1LPT = 1, 


OO OO 2 

-Pdm(fc) = [l + «lLPT(fe)^^] E -PsPT, 2 —loop(^) EE -fsPT,(i-i)- loop«i,lLPTA;^-^ (2.6) 

i=0 i=0 j=0 

The lowest orders (np to the 2-loop) are, 

p’dm(^) = Phik) + PsPT,l-loop(^) + aipLPT^^P’L(^) 

+ -PsPT, 2-loop (^) + aipLPT^^PbPT, 1-loop (fc) + a2,lLPTfc^pL(^) + • • • (2.7) 

At 1-loop SPT order the leading correction is P’L(fc)cii,iLPTfc^) same as in the standard EFT 
approach [2]. At 2-loop order we pick np two EFT terms in addition to the 2-loop SPT term, 
«i,iLPTfc^PsPT,i-loop(^) and a 2 ,iLPTfe^Pb(^)- If the radins of convergence for EFT is larger 
than for SPT then we may not need the latter term. Note that 2-loop SPT term vanishes 
at low k in 1-d relative to 1-loop SPT: this is no longer the case in 3-d, as discnssed in next 
section. 

We nse the coefficients derived above and SPT terms from [22] to plot the error of 
the different expansions relative to the fnll solntion. The resnlts are shown in fignre 1. 
Without any EFT parameters the ILPT (and 1-loop SPT) solntion has 1% accnracy only to 
k = 0.03h/Mpc. With 1 EFT parameter applied to fnll ILPT (i.e. infinite loop SPT) we 
find 1% accnracy to /c ~ 0.2h/Mpc (solid bine line in fignre 1), and with 3 EFT parameters 
to A: ~ 0.45h/Mpc (solid orange line in fignre 1). In contrast, 1-loop SPT applied to fnll 
«i,ilpt(A:) (i.e. infinite order EFT) is 1% accnrate to A: ~ 0.04h/Mpc, and 2-loop SPT to 
k ~ O.lh/Mpc. There is no improvement in adding additional EFT parameters np to this 
order. This confirms that the radins of convergence of transfer fnnction (EFT) expansion is 
larger than SPT radins of convergence. Beyond 2-loop the improvements in SPT loops are 
more modest: 5 loop SPT only extends the agreement by 10% in A:, mostly becanse 5 loop 
SPT does not improve mnch the agreement with ILPT. Overall the EFT gains in combination 
with SPT are most snccessfnl at very low k, where one expects one EFT term ccipLPT to be 
snfficient and SPT loop terms are rapidly converging to ILPT. 

There are several lessons of 1-d example worth emphasizing. As shown in [22] in 1-d case 
ILPT is exact solntion at the pertnrbative level since the eqnation for the displacement field 
is linear, bnt shell crossings invalidate this solntion: hence PT can never fnlly describe com¬ 
pletely the dynamics of dark matter. The SPT series can be resnmmed into ILPT solntion, 
so ILPT is snperior to SPT at any given loop order. One can define a concept of trans¬ 
fer fnnctions T^^Lprp(A:) that is defined as a ratio of dark matter to ILPT power spectrnm, 
that contains all the information on the power spectrnm effects beyond PT. One can expand 
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T^^lpt(A:) into a series of powers of k (at low k only even powers contribute) , with coefficients 
that act as EFT parameters. This EFT expansion resummed and applied to ILPT then gives 
the full dark matter power spectrum. There are thus two expansions, one in terms of EFT 
parameters multiplying even powers of k, and the second expansion is related to the loop 
orders of SPT, and the latter has a shorter radius of convergence. The same EFT expansion 
of the transfer function derived from ILPT also applies to each of the SPT loop terms: there 
is only one EFT expansion in 1-d. One can order the two expansions in terms of EFT+SPT 
order. At 1-loop order one hnds the usual EFT expression, with the hrst EFT term (A:^ term) 
multiplying Pl(A:), while at 2-loop order the hrst EFT term also multiplies 1-loop SPT term. 
Finally, we note that 1-loop SPT+EFT does not extend the range of 1-loop SPT signihcantly, 
only by 30%, while 2-loop SPT shows a considerable improvement, by a factor of 3 in scale. 
A third scale one can dehne is the scale where stochastic terms become important, and where 
ILPT no longer correlates well with dark matter. There is no obvious advantage in separating 
the terms into the part that correlates with ILPT and the part that does not, if one is only 
interested in the power spectrum. This changes if one also includes higher order correlations, 
as discussed below. Many of these features translate into 3-d as well, but there are additional 
complications in 3-d that are described in the next section. 

3 EFT expansions in 3-d 

We have seen than in 1-d the scale dependence of r^LPT(^) with even powers 

of k and that by expanding Pilpt into an SPT loop series there is a well dehned procedure 
that gives us an expansion in both EFT parameters and SPT loop order. As already shown in 
[17, 18, 23], in 3-d there is also very little stochasticity between 2LPT and dark matter at low 
k, about 10% in terms of total contribution to the overall nonlinear effect for k < 0.2h/Mpc. 
For ILPT the stochasticity is larger, and cannot be neglected relative to other nonlinear effects 
at any k. This is similar to the 1-d case and has the same origin: ILPT and 2LPT determine 
the positions of halo formation (where shell crossings occur) well, but LPT displacements do 
not stop there, but instead particles continue to stream and spread the dark matter by a 
distance that is approximately the same everywhere, independent of the halo mass. In terms 
of the EFT parameters one therefore expects at the lowest order a similar correction as in 1-d 
case, which is of order of Q:iLPT^^PL(^)) where cCi/px is several Mpc/h and the correction is 
positive relative to ILPT or 2LPT. In general, we can again dehne 

Pdm = PiLPT (l + aiLPT(fe)A;^) , (3.1) 

where for this paper we have i = 1,2, 3. We will also apply this to CLPTs model presented 
in [20], where the linear power spectrum is truncated at the nonlinear scale. In the following 
we will use P^m = P^lpt-^slpt using T 2 lpt from simulations, since 2LPT correlates very 
well with the dark matter and 1 — < 1% for k ~ 0.2h/Mpc (and much less than that for 

lower k [17]). This ensures that sampling variance cancels to a large extent, and since P 2 LPT 
can be computed analytically we thus have a measurement of Pdm without the sampling 
variance down to very low k. Our analytic 2LPT power spectrum, P 2 LPT) is based on the 
1-loop LPT calculations presented in [20]. As it was shown there, in 2LPT case, second-order 
displacement is used to compute the corresponding density power spectrum, and since the 
relation between the density and displacement is nonlinear, perturbative analysis relays on 
the cumulant expansion. Truncation of this expansion then implies that 2-loop and higher 
LPT terms that show up in higher cumulants are ignored. Nevertheless, this approximation 
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Figure 2. Running of a{k) for severals different models. On the left panel we show the running of 
the LPT models (in green) related to Eq. (3.1): ILPT (solid line), 2LPT (dashed line), and 3LPT 
(dot-dashed line). We also show the CLPTs model in Eq. (3.13) (purple solid line). On the same 
panel we show the running of the a’s related to the hybrid models from Eq. (3.15): Hyl (blue solid 
line), Hy2 (orange solid line), Hy3 (blue dashed line) and Hy4 (orange dashed line). On the right 
panel we show one loop (red dashed line) and two loop (orange dashed line) results for SPT EFT 
models, and also the IR resummed verisins of the same lines (solid red and orange lines). One loop 
(blue solid line) and two loop (blue dashed line) Hyl results also shown, as well as one loop results of 
LEFT [25]. 


was found to be 1% accurate against simulations for scales k < 0.14h/Mpc [20], and above 
this scale we switch from analytic model to the N-body simulation results. 

We performed simulations using wiggle and no wiggle realizations of the same initial 
conditions seeds, using fastPM code [24]. We construct the no wiggle linear power spectrum 
from the wiggle power spectrum using the method first introduced in in [20] and explained 
in detail in appendix A, which ensures that these two have the same ag, as well as velocity 
dispersion (t„. For this work, flat ACDM model is assumed with Om = 0.272, LIa = 0.728, 
rib/^m = 0.167, h = 0.704, rig = 0.967, ag = 0.81. The primordial density field is generated 
using the matter transfer function by CAMB. We ran several simulations of both 1.3Gpc/h 
and 2.6Gpc/h box size and 2048^ particles, each with a wiggle and no wiggle initial conditions. 

Next we consider the dependence of a parameter on scale k. In the case when a{k) 
shows any scale dependence we call it running of a{k) parameter. The motivation is to test 
these models and to determine which one exhibits the least running of a{k) over some range 
of scales. The absence of running improves the predictive value of a given model since only 
a constant parameter needs to be determined while in the case of running the full functional 
form is required (over the same range of scales). For example, in EFT framework, we expect 
a’s to be constant parameters so in this picture observed scale dependence at a given scale is 
the indication of importance of higher order corrections that have not been included. 

In figure 2 we present the running of aiLPT(^) as a function of k, analogous to 1-d 
analysis (top of figure 1). We only plot the results for BAG wiggle case, as the wiggle and no 
wiggle simulations give near identical results for the transfer functions. Note that aiLPT(fc) 
show no evidence of BAG features: most of the BAG information is thus in Pilpt already. We 
see considerable running of EFT parameters aiLPT(^) over the entire range of k: we note that 
even a small amount of running at high k can lead to large effects in Pdm(^)- The accuracy 
of these measurements depends on the accuracy of simulation data. Gn scales smaller than 
k ~ 0.08/i/Mpc simulation power spectrum can be taken as 1% accurate (see [24] for details). 
At low k, on the other hand, strong running is exhibited that is most likely caused by the 
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Figure 3. Scale dependence of the linear two point functions of displacement field, which contribute 
to the cumulant expansion, Eq. (3.6). We have split the contributions into the wiggle (left panel) and 
no-wiggle part (right panel). All results are shown at redshift z = 0.0. We see that the wiggle part 
has a most of the support at scale ~100 Mpc//i. 


inaccuracies in simulations related to the used speedup methods, (see [24] for details). 


3.1 BAO damping derivation 

Recently several resummation procedure have been suggested that include the effects of linear 
displacement field [4, 25, 26] on two point statistics of wiggles. Such resummations do not 
improve the reach of the perturbative expansions but they do correctly predict BAO damping. 
BAO effects show up in the power spectrum as 2% residual amplitude oscillations. In [26] 
a simple approximate procedure has been presented that captures the bulk of the damping 
effect on the BAO wiggles. Here we provide an alternative derivation. 

We start by dividing the initial power spectrum into the non-wiggle and wiggle part: 


— ^nw,L (A:) + AP«,L(fe) 


(3.2) 


where Pnw it the broad band part and AP„ is the residual containing only the wiggle part. 
This split, given above, is not unique which allows us to impose further constraints on both 
Pnwik) and AP^. We requite that the resulting nonlinear power spectrum obtained form to 
Pl,w and Pl.dw give us the same broad band power amplitudes on large and small scales. 
This we can achieved by requiring that erg and ay be the same for the Pl,w and PL.nw power 
spectrum, i.e. 

^8 = J W‘^{kR)PL,^{k) = 0 and (^v = \j (3-3) 

where ITr(A:) = 3 [sm{kR)/(kR)^ — cos{kR)/(kR)"^] and i? = 8Mpc//i. Detailed construction 
of the Pnw using the constraints above is given in Appendix A. 

In Lagrangian formalism the power spectrum can be written as [20, 25] 


{2Trf6^{k) + P{k) = I d-^qe 


-*q-k /g-ik-A\ ^ 


= d^qe 


—iqk 


exp 


-\kikjAij{c\) + ... 


, (3.4) 


where dots “... ” represent the three and higher point cumulant contributions, and we have 

Aij(q) = (AiAj)^ = X{q)Sfj +Y{q)qiqj (3.5) 
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where the two scalar functions are given by 


poo 

X{<1) = / 

poo 

yiQ) = / 

Jo 




TT 


ji(feg) 

kq 


dh 

^PmMkq), 


IT 


(3.6) 


and Pip{k) is the diagonal part of the displacement power spectrnm [20]. In linear approx¬ 
imation we have Pijj{k) —)• Pi^{k). Since X and Y are linearly related to we can simply 
separate 


— -^ij,nwiq) P ( 3 -'^) 

Keeping only the linear part of Aij exponentiated and expanding the rest we can rewrite the 
power spectrnm 


PAk) 


j d?q e 2 K^i^L,nw(‘i) _|_ “higher no — wiggle terms’’^ 

-I- J d^q ^ _ ^kikjAA^^^{q) + “higher wiggle terms’’^ 

Pn^k) + J d\ _ hikjAAlJq) + ... ) (3.8) 


where “... ” now represent the terms of one loop order and higher, and Pnw is the no-wiggle 
nonlinear power spectrnm given by the same expression as in Eq. (3.4), with the initial power 
spectrnm Pn.nw Expanding the angnlar part of the ^Lnw keeping the monopole part 
exponentiated we have 




- -A:i%A^Jf^^(q) + ... 


= f d^q 2*^^('’*^L.nw+3>L,nw) 


= I d\ 


- “ Ihkjimj - l/3h,^)yL,nw) + 


= e 




^kikjXA^^^P ^k kikjAA-^^^P2{k')yL,nwP 
APh,w + “higher order wiggle terms” .. 


(3.9) 


where V 2 {k) is the second Legendre polynomial. We eventnally neglect the term propor¬ 
tional to the V 2 {lP) angnlar dependence which, at the leading order, leaves only hrst term 
in the sqnared brackets. We have also introdnced averaged qnantity T? = ^^^AL,nw(9) + 
^lL,nw(9)))^ — 28.0(Mpc//i)^ where ((•)) represents the averaging over q range where snpport 
of AAJ^^ is prominent, i.e. q ~ 110Mpc//i. In Eq. (3.9) above we have nsed the fact that 
in the region where has a non-negligible snpport (aronnd the scale of ~ 110Mpc//i), 
varies slowly. This can be seen in hgnre 3 where X and Y components of both wiggle 
and non-wiggle part are shown. This approximation also known as Laplace’s approximate 
integration method. We note that fnlly nonlinear X and Y are nearly identical to the linear 
aronnd q ~ 110Mpc//i [27]. It has been argned [4] that this IR resnmmation does not affect 
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the broadband SPT or EFT terms, so we will leave the broadband part nnchanged. Thns onr 
fnll IR-SPT model at one loop gives 

-Pdm(^) ~ Pnw,L(^) T Pnw,SPT,l—loop(^) T ClSPT,l—1 oop,Ir(^)^ Pnw,L(^) (3.10) 

+ e~^ ^ AP„^sPT, 1-loop(fc) + (l + (“SPT,1-loop,IR + AP„^l(^)) ■ 

where we have introdnced one-loop wiggle only power spectra 

^-^WjSPT,!—loop (^) -fw,SPT,1—loop (^) -fnw,SPT,l—loop (^) ■ (^'H) 

Similar procednre can also be straightforwardly applied to the hybrid and 2-loop SPT resnlts 
as discnssed farther below (see also [18]). While we have focnsed on BAO wiggle here, the 
derivation applies to any wiggle localized in q. In general, the lower valne of q the lower the 
valne of damping distance S. 

3.2 Power spectrum results 

We observe that the running for ILPT and 2LPT is larger than in the 1-d case. The main 
difference between 3-d and 1-d analysis is that for the latter the full displacement solution is 
given by the first order ILPT (Zeldovich) displacement ^i, while in 3-d one has an infinite 
series of displacements, 'J' = P+T'i -|- D'^^2 + -Djj.T'a -|- .... At 1-loop SPT order in P{k), we 
have contributions up to 3LPT (^'s). But this term is not included in P 2 LPT) which therefore 
does not contain the full 1-loop SPT. Similarly, ILPT does not include both ^^2 and ’^' 3 . 
If, as expected, the correct dark matter solution contains the full SPT 1-loop terms, then 
the running of aiLPT(fc) and a 2 LPT(^) reflects the k dependence of the difference between 
ILPT and 2LPT at 1-loop level versus the full 1-loop SPT. In contrast, 3LPT contains all 
1-loop SPT terms and we expect the running of q; 3 lpt(^) to have less k dependence at low k. 
However, 3LPT contains additional terms beyond 1-loop SPT that are quite large and largely 
spurious [20], so we do not expect 3LPT EFT parameter to agree with the corresponding 
SPT 1-loop EFT parameter defined below. For example, at low k 3LPT has an additional 
zero lag contribution k'^Pi^[k)af^, which is a 2-loop SPT term, and is quite large. At higher k 
there are additional 2-loop SPT terms in 3LPT, which cause q: 3 lpt(^) to be running with k. 

To gain more insight into the running of these parameters we can perform the expansion 
of 1-3LPT at low k, which has been shown to be accurate for k < 0.12h/Mpc [20]. We have 

PsPT,i-ioop(A:) = (1 - k^al)PUk) + ^Qsik) + ^Qi{k) + ^i?i(fe) + ^ (Q2{k) + 2 R 2 {k)) , 
PiLPT,i-\oop{k) = (1 - k‘^crl)PL{k) + -Qsik), 

P2LPT,i-ioop(fe) = (1 - k^^l)PUk) + ^Qsik) + ^(^Q^{k) + 2R2ik)) + ^Qi(fe), 

-P3LPT,i-ioop(A:) = PsPT,i-ioop(fe) - k'^PL{k)alioop, (3.12) 

where we have introduced several LPT terms defined in e.g. [28]. However, the expression 
for p 3 LPT,i-ioop is not necessarily valid at low k. The reason is that there is a large zero 
lag correlation of 22 and especially of 13 displacements. At the next order in 3LPT we have 
-PsLPT = P’sPT,i-loop(fe) - k‘^PL{k)aly^^p, where is the sum of zero lag correlations of 

22 and 13 LPT displacements. At 2 ; = 0 the linear theory value of r\j 36 (Mpc//i)^, and 
'^iloop ~ (Mpc//i)^, so the 2-loop SPT term in 3LPT is quite large even at very low k. It 
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is almost entirely spurious [20] (see also [18] for recent results). In [20] we introduced CLPTs 
model that attempts to cure this problem by truncating high k contributions. We can define 
a model related to the CLPTs in a similar way as in Eq. (3.1) 

^’dm = .PcLPTs (l + acLPTs(fc)A:^) , (3.13) 

Results for this model are also shown in figure 2. We observe that EFT term is close to 0 
around k ~ O.lh/Mpc, but it still has a strong k dependence. 

We see that the 2LPT solution does not contain all the terms at the 1 loop SPT level. 
Why is this not important for its correlation with full dark matter, i.e. why can we drop 1-loop 
SPT term in 2LPT and still have near perfect correlation with the final dark matter? This 
is because the missing term is dominated by Ri{k) ~ aji^{k)k‘^Pi,{k). Since we are allowing 
a fully general a 2 LPT(fc)) it can include this term running with k. Since aii^{k) is strongly 
changing with k even for k < O.lh/Mpc, it leads to a strong running of Q; 2 LPT(fc)- Figure 
2 shows there is also considerable running of 3LPT, even though it contains all 1-loop SPT 
terms. This suggests that 2-loop and higher order contributions cannot be neglected except 
at very low k, as also shown in [20]. These higher order terms may however be spurious. 

For the SPT perturbative expansion we may expect EFT parameter to be running less 
at low k, just as in 1-d case, at least in the range where 1 loop SPT is expected to be valid. 
We can define 1 loop SPT EFT parameter asPT,i-ioop as [2] 

Pdmik) = D\PL{k) [l -I- Q;sPT,l-loop(fc)/c^] -I- Ll+P’sPT,l-loop, (3-14) 

and derive the running of asPT,i-loop(^) from it. The leading EFT terms scale with Pi^{k)^ and 
hence contain BAO wiggles. It has been shown that this SPT EFT version is a poor model in 
the context of BAO wiggles [2] . At higher k the higher order effects of long wavelength modes 
leads to suppression of these BAO wiggles [4]. This resummation procedure is automatic in 
LPT schemes, which are therfore superior for BAO damping. 

Figure 2 suggests that aspT,i-loop,lR(^) is nearly constant for k < O.lh/Mpc, with a 
value around -2(Mpc/h)^. This value increases to about -4(Mpc/h)^ for k > O.lh/Mpc. This 
change suggests that the value of asPT,i-ioop,iR for ^ > O.lh/Mpc is no longer the low k EFT 
parameter at the lowest order. We will test this model further using BAO oscillations in next 
section. 

In 1-d we have identified LPT approach as clearly superior to SPT approach because 
of the automatic resummation of all SPT terms. Doing SPT loops has no advantages over 
ILPT, since only in the infinite loop limit it converges to ILPT, and the convergence is slow. 
In 3-d this is no longer the case, because higher order LPT displacements become less and less 
reliable, and their resummation is not necessarily a good thing. However, doing resummation 
on ILPT only is still likely to be useful, since it is dominated by modes in the linear regime, 
hence reliably computed by ILPT resummation. In 3-d, 1-loop SPT identifies properly all 
1-loop terms and is not contaminated by 2 and higher loop orders, and as a consequence we 
expect it to give a constant EFT parameter over the range of its validity and over the range 
of constant a EFT parameter validity. Moreover, we expect the EFT parameter relative to 
1-loop SPT to be small [2], and indeed the values in figure 2 suggest it is a lot smaller than 
the corresponding LPT versions. 

We thus want a scheme where we resum only ILPT (or ILPT and 2LPT), and add the 
remaining SPT 1-loop terms, and finally add or multiply EFT terms to it. We can define 


11 


several versions of this proposal, 

Hyl: Pdm(^) = -PlLPt(^) (l + aiLPT,l-loop(^)^^) + (^-PsPT,l-loop(^) — -PlLPT,l-loop(^)^^^ 

Hy2: Pdmik) = -P2LPt(^) (l + a2LPT,l-loop(^)^^) + (^-PsPT,l-loop(^) - -f2LPT, 1-loop 

Hy3: Pdm{k) = -PiLPt(^) + (-PsPT,l-loop(^) - -PlLPT,l-loop(fc) + PLO:SPT,lLPT,l-loop{k)k'^] . 

(3.15) 

All these models can readily be derived from a full resummed form of the power spectrum 
given in Eq. (3.4). The difference is only on which terms we keep exponentiated and which we 
expand. In the case of Hyl model ILPT displacement (Zel’dovich part) remains resummed 
while rest of the one-loop contribution is expanded. Similarly for Hy2 model we leave 2LPT 
displacement contributions resummed while we expand the one-loop residual. The difference 
of Hyl and Hy3 models is in the way contributions proportional to a is resummed; in Hyl case, 
it is proportional to Zel’dovich contribution and in Hy3 case to the linear power spectrum. 
These results are also shown in hgure 2 and show considerable running of EFT parameters. 
Models Hyl and Hy3 are almost equal, and we will show below they give very similar results. 
The same is also true for recently developed LEFT [25], which as well shows substantial 
running of a{k) at all k. 

In 1-d we have seen that SPT expansion have a smaller radius of convergence than the 
EFT expansion, so to improve the model we could repeat the EFT procedure on 2-loop SPT 
even with just 1 EFT parameter. We can assume that 1-loop SPT EFT parameter remains 
unchanged. Unlike the 1-d case, we must also absorb with the EFT parameter the very large 
zero lag value of which contributes at 2 loop order a term —A:^crgprp 2 -\oop^Lik), which 

dominates at low k, but which is almost entirely spurious [18, 20]. However, we would have 
to also absorb its 2-loop counter terms. It is clear that doing 2-loop SPT is considerably more 
complicated than in 1-d. If we assume low k 2-loop SPT is completely spurious at low k this 
is equivalent to the requirement that 2-loop SPT vanishes at low k, just as in 1-d case, 

Pdm{k) = (l + aSPT,l-loop(^)fe^)-PL(^) + (l + «SPT,l-loop(^)^^)HspT,l-loop(^) 

+ (-PsPT, 2-loop (^) — ^^<7sPT,21oop-PL(fc))- (3.16) 

Result for the running of a for this model is shown in hgure 2. More generally, if we split 
UgpT 2 ioop ^ spurious part and a real part then one hnds one needs a different EFT term 

multiplying Pp and PsPT,i-loop, 

Hdm(^) = (l + asPT,2-loop(^)fc^)RL(fc)-|-(l-|-asPT,l-loop(^)fc^)RsPT,l-loop(^)+RsPT,2-loop(^)- 

(3.17) 

In 1-d we have seen that at 2-loop level one can assume Q!SPT, 2 -ioop(fe) = «spT,i-ioop(fe) and 
that one does not need to add its /c^ dependence, but in principle one could add that as 
another parameter. This has been explored in [6], but we do not pursue it further here. 

3.3 Higher order correlations 

So far our approach has been to quantify the ignorance of PT by the transfer function, which 
can be translated into an effective EFT parameter that is htted to simulations. Over a 
narrow range of k this is guaranteed to give correct answer, and hence the method is simply 
parameterizing the ignorance with essentially no physics. However, the concept of transfer 
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Figure 4. BAO wiggles, i.e. ratio of the wiggle and non-wiggle power spectrum, is shown for four 
different models: Acdm (top left panel), Monodromy model (see [29, 30]) (top right panel) and two 
other models labeled V 3 and V 4 with additional wiggles relative to Acdm (bottom panels). Linear 
theory results (blue lines) evolve due to nonlinearities and yield results given by N-body simulations 
(black points). Wiggle damping for all these models is well described by the ILPT (Zel’dovich) 
model (green dashed line). Note that all the initial wiggles are highly dampened at lower scales, 
k < 0.5/i/Mpc. All results are shown at redshift z = 0.0. 


functions in the absence of stochasticity can be more useful than that: it allows one to derive 
the higher order correlations in the regime where r 2 LPT ~ 1- For example, 2LPT already 
contains all of the terms that determine density perturbation at 2nd order, 

52LPT(k) = 6Uk) + j d3kid3k2F2(ki,k2)<5^(k- ki - k2)<5L(ki)5z,(k2), (3.18) 

where F 2 (ki,k 2 ) is the SPT 2nd order kernel (see e.g. [8]), 

«(k,.k.) = F + i (I + I) k. . fe + 1 (3(k. . fey - 1 ) , (3.19) 

where k* = kj/fep Bispectrum is defined as 

(27r)35^(ki + k2 + k3)i?(ki, k2, ks) = {6{ki)6{k2)5{k^)), (3.20) 

and hence, at the tree level, one can write 

5(ki,k2,k3) = T{ki)T{k2)T{k3)(^F2{ki,k2)PL{ki)PL{k2) + 2 cycl. perm.^ (3.21) 


At the lowest order we can write T{k) = (1 -|- Q; 2 LPT(^)fc^)^'^^ ~ 1 -|- a 2 LPT{k)k‘^/2 and thus 
the above expression becomes 


-B(ki, k2, ks) = ^1 -b Q:2 LPt(^i)^i/2 + a2LPT(^2)^2/2 + a2LPT(^3)^3/2^ 

X F2{ki,k2)Ph{ki)Ph{k2) -b 2 cycl. perm. 
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To this one needs to add 1-loop bispectrum from 2LPT terms, which are at the same order 
as the EFT corrections of tree level bispectrum. This requires expanding <52 lpt to 4th order 
in (5 l. The 2LPT 1-loop bispectrum contains some, but not all of the terms of the full SPT 
bispectrum at the 1-loop order. Some of the missing terms may cancel out the scale depen¬ 
dence of Q! 2 LPt(^) in the EFT corrected tree level bispectrum above, but at the operational 
level this is not relevant: the expression above should include all of the terms at 1-loop level, 
including the EFT terms and should be valid as long as r 2 LPT = 1- As mentioned above 
stochasticity in 2LPT has about 10% contribution to the nonlinear terms for k < 0.2h/Mpc, 
and rapidly grows above that, so the expressions above should only be valid at this level. 
Previous analyses [31, 32] have argued that there are 4 different EFT parameters (although 
only 1 matters for improving the fits), while in the expression above there is only the transfer 
function, which is the same that also enters into the power spectrum calculations. While 
this transfer function, when expressed in EFT terms q: 2 LPt(^) runs with k and cannot be 
approximated as a constant, one can derive this running from the power spectrum. As long 
as we work up to the same k, and stochasticity can be neglected, there are no additional 
EFT parameters for bispectrum relative to the power spectrum. At the leading order any 
EFT corrections can only be a multiplicative factor times F 2 kernel, analogous to the power 
spectrum situation, where at the leading order one can only have EFT corrections multiplying 
the linear power spectrum PL{k). This may help explain why the additional parameters in 
EFT fits of bispectrum were found not to be needed [31, 32]. Same concepts can be applied 
to higher order correlators (trispectrum etc.) as well. We do not pursue this approach further 
here, but it would be interesting to see how it compares against the standard EFT bispectrum 
calculations [31, 32]. 

4 BAO residuals in EFT schemes 

Armed with these various EFT expansions we can test them on BAO wiggles. Let us first 
focus on the BAO residuals in standard Acdm model. The main idea is that EFT terms scale 
with Pl(^)) and hence contain BAO wiggles, and the better the scheme the better it should 
be able to explain the BAO wiggles. We compute the power spectrum for both wiggle Pw{k) 
and no-wiggle P^^{k) simulations and take the ratio of the two. This is shown in figure 4 top 
left panel, together with linear theory and Zeldovich approximation (ILPT). It is clear that 
ILPT does a good job in describing the BAO damping. In the same figure we explore also 
the wiggles of some other (non-BAO/Acdm) models, which are discussed further below. 

Next, we look at the small wiggle residuals. In order to do this subtract out the ILPT 
part and then compare the models and N-body simulations. This is shown in figure 5. We 
see that there is a residual BAO, even if it is small, about a factor of 10 smaller than the 
original BAO. It is these residuals that we wish to test against the EFT expansions. Our 
model for BAO residuals is simply given by P’dm,w(^) and Psm nw(fc) versions of the EFT 
models presented in previous section. EFT of ILPT ^ is a decent model and predicts some of 
the residuals, but 2LPT and 3LPT (all defined by Eq. (3.1)) are a lot better. For ILPT the 
stochastic term is not negligible [17], so we expect a worse agreement if we apply Q:iLPT(fc) 
to model BAO residuals. In the same family we can also add the CLPTs model (see [20] 
) for which EFT extension is defined in a same way as it was for iLPT models, Eq. (3.1). 
Performance of CLPTs on the wiggle residuals is of the same level of accuracy as 3LPT. Next, 
we look at the hybrid models defined in Eq (3.15). The difference of these models are in 

^We note that our ILPT EFT model is equivalent to the ZEFT model introduced in [25]. 
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Figure 5. Residual wiggles, relative to the ILPT (Zel’dovich). In top two lines iLPT models, as 
well as CLPTs model (see [20]) are shown using definitions in Eq. (3.1). In panels in lines three and 
four we show residuals of hybrid models defined in Eq. (3.15), as well as LEFT theory developed 
recently in [25] . In a bottom line we show SPT-EFT one loop and two loop models (IR resummation 
included), given by Eq. (3.11) and (3.16) respectively. All results are shown at redshift z = 0.0. 


resummation of up to ILPT terms (for Hyl and Hy3 models) or up to 2LPT terms (for Hy2). 
These models can be considered to be at 1-loop level (up to the resumation). We see that 
ILPT resummation models (Hyl and Hy3) work better on the wiggle residuals then 2LPT 
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(Hy2). Indeed residuals of Hyl and Hy3 are amongst the best models considered and their 
performance is matched only by two-loop SPT model. LEFT [25] performs at the similar level 
of accuracy as IR resummed one-loop SPT. For one-loop SPT we also find a decent fit for 
IR versions, and not very good fit for non-IR versions (as can be already seen from figure 2). 
Going to two-loops IR resummed SPT improves the wiggle residuals further and gives a very 
good agreement. Thus our results suggest that these residuals can be modelled as a smooth 
function multiplying PT power spectrum. These results suggests that the broadband analysis 
extracting transfer functions relative to LPT or SPT is also able to reproduce the BAO 
wiggles, so the picture is consistent, and in all cases introducing EFT parameters improves 
the agreement on the wiggle part of the power spectrum. 

4.1 Wiggles in primordial power spectrum 

Here we focus on the wiggles and wiggle-like features in the power spectrum beyond the 
Acdm. These could be imprinted by the physical processes during inflation, for example. 
It is interesting to see how the models discussed above perform in predicting such features 
given that a(k) has been determined from the broad band spectra. For this purpose we 
construct three new wiggle models; first is the monodromy-like model (see [29] and [30] for 
the parameteization we have adopted), which has oscillations in Ink, and the other two 
models (labeled V 3 and V 4 ) are using the same BAO wiggle power spectrum with the boosted 
amplitude and shifted scale dependence. These additional wiggle power spectra are then 
added to the Acdm power spectrum. In figure 4 we show the ratio of the total (Acdm plus the 
additional wiggles) wiggle to non-wiggle power spectrum. Linear theory shows the features 
of the initial power spectrum which are then subject to the nonlinear evolution. We can see 
that all the initial wiggles are highly damped at high k, k > 0.5/i/Mpc, in fact, wiggles of V 3 
model are completely washed our at redshift 2 ; = 0 . 0 . 

We then again look at the small wiggle residuals for these models, where we subtract 
out the ILPT part and then compare the models and N-body simulations. This is shown in 
figure 6 for monodromy-like model and in figure 7 for V 3 model (as V 4 residuals are equivalent 
to the Acdm at 2 ; = 0.0 we do not show them). In addition to the ILPT residuals we show 
the 2LPT, EFT-SPT and Hyl model residuals. We see that as in the earlier case EFT-SPT 
and Hyl perform very well and reproduce the wiggle shape to high accuracy. 

5 Interpreting the dark matter power spectrum 

In this paper we have addressed the EFT modeling of power spectrum, which introduces EFT 
parameters that parametrize the ignorance of PT and at low k scale as multiplied by a 
low order PT. These have to be supplemented by stochastic terms at higher k. In [16, 33] 
a halo inspired model for the dark matter clustering was advocated, where one takes ILPT 
(Zel’dovich) for the 2-halo term, and adds to it an effective 1-halo term that is very localised 
(to a few Mpc/h). It has been shown that one can build a model with high accuracy using 
this kind of ansatz. The 1-halo term is not necessarily defined as the mass within the virial 
radius: there is nonlinear clustering outside the virial radius that adds to this effective 1 -halo 
term. In this section we look at the relation between the two approaches. 

Both 1-d and in 3-d LPT simulations suggest that ILPT (and higher order LPT) tem¬ 
porarily create dark matter halos in the regions of orbit crossings, but in LPT particles 
continue to stream through. In full dark matter simulations particles stick together after or¬ 
bit crossing, locked inside the high density regions called halos. In LPT the particles continue 
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Figure 6. Residual wiggles for the monodromy model shown in 4, relative to the ILPT (Zel’dovich). 
We compare the performance of transfer function extended of ILPT (green dashed line) as well as 
2LPT (green solid line) models (given by Eq. (3.1)). In addition we show the performance of one 
loop (red solid line) and two loop (red dashed line) SPT models as well as one loop Hyl model (see 
Eg. (3.15)). Note that the same values of a{k) are used as in A model. Results are shown for redshift 
z = 0.0. 


to stream out of these high density regions, effectively smearing these halos of to a few Mpc/h. 
As a result the power spectrum in any LPT is below the linear power spectrum (at z = 0). 
However, the seeds of halo formation have been imprinted already at the LPT level, and show 
up in the transfer functions, which essentially account for the artihcial smearing of the high 
density regions by a smoothing radius that is independent of the position and hence can be 
represented as a k dependent transfer function. Multiplying LPT power spectrum with the 
square of the transfer function T{k) reverses this smearing. The halos, at least massive ones, 
exist in LPT, but their mass is spread out to such a large radius that they are not visible 
except for the largest ones. Reversing this with the transfer function brings the total halo 
mass to a smaller radius, and produces a better dehned 1-halo term. These halos are still 
too diffuse relative to the N-body simulations, but contain all of the mass that has collapsed. 
As a consequence this term gives the correct 1-halo amplitude at low k {k ^ 0.2h/Mpc for 
z = 0). 

The transfer function T{k) only goes so far in improving the relation between LPT and 
full dark matter N-body simulation. On smaller scales the dependence of the halo prohle 
on the halo mass becomes important: halos have internal prohles, which depend on the halo 
mass, with larger halos having a larger extent: the relation can no longer be described in 
terms of a spatially independent transfer function. This is the stochasticity (sometimes called 
mode coupling) term Pj{k). To improve the relation between LPT and N-body simulations 
one needs to further compress the halo prohles of the halo blobs as dehned by the transfer 
function corrected LPT to the true halo prohles. The typical scale for this corresponds to the 
typical halo radius squared, averaged over the halo mass, which is of the order of IMpc/h, 
a smaller scale than the LPT smearing scale inside the transfer function. This process only 
redistributes the halo mass while conserving the total collapsed mass, so the process only 
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Figure 7. Same as figure 6, for the model V 4 shown in figure 4. 


becomes important at higher k and does not affect clustering at k < 0.2h/Mpc. 

Note that both the transfer function correction term and the stochasticity term essen¬ 
tially try to do the same thing: correct the overly diffuse nature of halos in LPT to a more 
compact form found in N-body simulations. The only difference is that the transfer function 
part accounts for the halo mass (or spatial position) independent part, and dominates at low 
k, while stochasticity term depends on halo mass, and hence spatial position, and dominates 
at high k. They are also both determined by terms beyond PT, and become essentially free 
functions at high k. For the power spectrum there is therefore no obvious reason for splitting 
effects beyond PT into the two terms, and indeed in this paper we also defined T(k) as a 
combined effect. Similarly, in [33] both terms are put together into a single effective 1-halo 
PBB{k) term, 

-Pdm = -PlLPT + {TiBBT ~ l)-f’lLPT = AlPT + (T'iLPT “ l)-PlLPT + -Pj = -PlLPT + Pfes, (5.1) 

where PBB{k) is the non-perturbative part, that puts together all of the terms that cannot be 
computed in ILPT, which in [16] was modeled as Pbb(^) = P{k)-A-o{t — P^j^ ^k'^P---)- 

Here Aq is the amplitude of effective 1-halo term and Rih, 2 i are connected to the various 
moments of halo radius averaged over the halo profile and halo mass function, typically of 
the order of IMpc/h. For consistency of the halo model one needs to compensate the 1- 
halo term at very low k, and F{k) is the compensation term, which in [33] was modeled 
as F{k) = 1 — 1/(1 -|- k'^R^), where R 25Mpc/h. This term is not associated with any 
physical scale, and is instead related to the ILPT streaming error encoded in the leading 
EFT parameter (figure 2). 

Since this compensation term is determined at low k, where stochasticity for 2LPT can 
be neglected, one can derive it from 2LPT transfer function. This is less true for ILPT, so the 
difference between ILPT and 2LPT can be viewed as the source of stochasticity for ILPT. The 
effect of compensation is typically below 1%, so quite small. Since we have measured this term 
without sampling variance from 2LPT we can look at its scale dependence at low k, which is 
usually sampling variance dominated. In figure 8 we plot this term against the F{k) model 


- 18 - 






































of [33]. We see that this term is only approximately modeled by the F{k) ansatz above at 
very low k, although in practical terms this does not make much difference. In particular, the 
small EFT generated BAO effects discussed in previous section are not captured by PBB^k) 
ansatz (which was also noted in [33]). We note however that our results at low k are only 
approximate, both because of stochasticity of 2LPT (at 10% level) and because our analytic 
.P 2 LPt(^) does not contain terms beyond 1-loop. The EFT term is similar to a 1-halo term, 
meaning that it is flat in k, around k ~ 0.1 — 0.2h/Mpc. This is because the linear power 
spectrum scales as k~‘^ in that range. This is where one can merge the two approaches. In 
figure 8 we show an attempt where we do this, defining maximum in Pdm ~ .PiLPT in EFT 
model to switch from EFT approach to halo model approach, which happens roughly at 
k ~ 0.2h/Mpc. EFT model of Eq. (3.11) uses 1 EFT parameter. 

In [33] both Aq and F{k) were extracted from SPT at low k (0.02h/Mpc < k < 
0.03h/Mpc), under the assumption that the EFT correction is small. In figure 8 we see 
that this correction is indeed very small there, but determining both parameters from a nar¬ 
row range of k induces degeneracies between them. A more robust approach would be to 
determine AQF{k) using simulations at higher k, or using EFT approach, where EFT param¬ 
eters are fitted to simulations. However, since the effect of F{k) is less than 1%, in practice 
one is simply trading one free parameter, Aq, for another free parameter, that of EFT model. 
It may be possible to connect EFT and halo model approaches in the regime of overlap where 
both are applicable. 


6 Conclusions 

The goal of the present analysis is to clarify some connections between the different PT 
schemes (LPT and SPT), the corresponding full dark matter solutions, and the halo formation 
process. In 1-D, ILPT is exact at the PT level, corresponds to the SPT at infinite loop order, 
but does not give exact solution because of sheet crossings (“halo” formation) that are beyond 
PT. One can parameterize the PT ignorance by a transfer function squared, defined as a ratio 
between true power spectrum and PT power spectrum, and one can expand this ignorance 
into a Taylor series which starts as 1 -|- ak"^, which corresponds to the EFT expansion at 
lower PT orders. At higher orders (beyond the two loop power spectrum) this equivalence 
does not hold any longer, and one would in principle need different transfer functions for 
different PT orders of the overdensity field. In addition, one can decompose the ignorance 
of PT into a part that is correlated with the true density, and a part that is not, called 
stochasticity. The latter dominates on small scales, and correlated part dominates on large 
scales. The advantage of this split is that the correlated part can be used to predict higher 
order correlations of dark matter from the corresponding higher order correlations of the PT 
field, with no extra parameters as long as stochasticity can be ignored. In 1-d this EFT 
expansion has a larger radius of convergence than the SPT series, and knowing the two allows 
one to order the combined SPT+EFT expansion. In practice, we find 1 parameter EFT 
suffices for the range of 2-loop SPT. The improvements at 1-loop and 2-loop are impressive 
at low k, but get progressively harder as we push to higher k where SPT series is very slowly 
convergent. Applying the EFT expansion to ILPT, which corresponds to infinite loop SPT, 
is always better than the SPT+EFT expansion at the same order. There is only one EFT 
expansion, and only one set of EFT parameters, parameterizing the ignorance that can be 
applied to both SPT and ILPT. 
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Figure 8. Difference of nonlinear dark matter power spectrum Pdm and Pilpt is shown for various 
models. We show N-body simulation results (as black points) and compare them to the transfer 
function approaches (red solid line). For comparison we show the HZPT model [33] (blue solid line) 
and the two-loop EFT results from Eq. (3.16) (orange solid line). For two-loop EFT results we added 
the (orange) band showing the results for aspT,i-ioop = —2.2 ± 0.4 [h/Mpc]^. Results are shown for 
redshift z = 0.0. 


In 3-d the situation is more complicated due to the lack of an exact PT solution. One 
has to compute higher order displacement fields, which can be computed in LPT, but their 
corresponding higher order LPT contributions become less and less reliable, and analytic 
expressions for the density field cannot be readily computed. It becomes more reliable to 
keep all the terms at a given PT order, which corresponds to f-loop SPT. For each LPT and 
SPT scheme there is a corresponding PT ignorance expansion in terms of EFT parameters 
that can be defined, and in general there is considerable running of EFT parameters as a 
function of k. For some schemes, such as 1-loop and 2-loop SPT, we expect there is little or 
no running at low k, and we indeed find that up to fe = O.lh/Mpc for 1-loop and k = 0.2h/Mpc 
for 2-loop. We also explored hybrid PT approaches, combining the full power spectrum of 
ILPT with the 1-loop or 2-loop SPT terms missing from ILPT, but this does not extend 
the constancy of EFT parameter to a higher k (and similarly for LEFT [25]). All wiggle 
models include damping of the wiggles by IR modes, and we provide a simple derivation of 
the damping, including specific application to BAO. 

Finding a single EFT parameter is certainly desirable from the point of view of having 
a small number of free parameters, but does not address the issue of correctness of the EFT 
approach. To address that we produced N-body simulations with and without wiggles in the 
power spectrum to test these different EFT approaches. We ask the question whether the 
addition of EFT parameter improves the BAO wiggle description over the corresponding PT 
case without EFT. In all cases we find that the EFT term is able to reproduce BAO wiggles 
better than the corresponding PT models, supporting the basic idea of EFT. However, in 
many of the models PT alone does quite a good job in describing the BAO, and the residuals 
that EFT terms improve are very small. For k > 0.2h/Mpc (at z = 0) the difficulties of 
having reliable LPT or SPT higher order calculations, the running of EFT parameters, and 
the stochasticity all contribute to diminishing returns, and it appears difficult to extend EFT 
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approach meaningfully beyond this k without iutroduciug additioual free parameters. We 
also discuss the couuectiou of EFT to the halo model [33], arguiug that there is a regime 
arouud k ~ 0.2h/Mpc where both descriptious are valid, allowiug oue to couuect the two. 

The coucept of the trausfer fuuctiou describiug the uouliuear effects beyoud PT without 
geueratiug stochasticity ou large scales, is a useful tool that eusures there is ouly oue such 
fuuctiou that cau be applied to correlatious at all orders. This fuuctiou must scale as k"^ at 
low k, where it is determiued by a siugle uumber that cau describe mauy differeut statistics. 
However, there is a fuudameutal difficulty iu applyiug this coucept to LSS, iu that it works best 
at the lowest k. But iu LSS these scales are of little iuterest to be modeled precisely beyoud 
the liuear theory, because uouliuear effects are small, aud the sampliug variauce errors, which 
scale as the iuverse square root of the uumber of modes, are large. The uouliuear effects 
are a few percent for k < O.lh/Mpc aud we are unlikely to reach this level of precision 
observationally in the near future, and for k > 0.2h/Mpc the EFT modeling with a single 
parameter breaks down at a sub-percent level, leaving 0.1h/Mpc< k <0.2h/Mpc as the range 
where these models can have their main applicability. This is also the range where primordial 
wiggles still have some signatures in the nonlinear power spectrum. These statements all 
apply to z = 0, and for higher redshifts the corresponding scale is shifted to higher k. 
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A Construction of the no-wiggle power spectra 

Here we summarize the methods of smoothing of the wiggles produced by baryon acoustic 
oscillations (BAG). We look at two filtering methods. The first employes the spherical Gaus¬ 
sian Liters on spherical 3D power spectra and the second employs a ID Gaussian Liter on the 
logarithmic scale. The latter turns out to be a superior method, recovering correctly small 
and large scale limits. Finally, we make use of the basis splines (B-spline) interpolation in 
order to achieve BAG smoothing, allowing us to impose additional constraints on the resul¬ 
tant power spectrum. This procedure enables us to construct no-wiggle spectra that have the 
same dispersions (t| and as the wiggle spectra. This is useful since we would want that our 
wiggle and no-wiggle power spectra exhibit same broadband nonlinear evolution. Additional 
advantages of the spline-based method are computational efficiency and automatization. 

We Lrst explore the smoothing of the power spectra using Gaussian Lltering. First we 
note that it is useful to divide the initial power spectra by some approximate (wiggle free) 
curve in order to reduce the amplitude range of the broad band power. After smoothing is 
preformed, we retrieve the total no-wiggle power spectrum by multiplying the result back by 
this approximate curve. Final no-wiggle power spectrum is then given as 

Prm{k) = Happrox(^)-^[-P(A:)/T’approx(fe)] , (A.l) 

^We use the approximation result given in [34], but alternative methods like BBKS [35] would work 
equally well. Alternatively, one can use either some power law form fe", with n oc 1.5, or actual output of the 
Boltzmann codes with no baryons. 


- 21 




Figure 9. Linear power spectrum is shown in dotted back line. Solid blue line is the smoothing result 
obtained by applying the ID Gaussian filter, and solid purple line is obtained applying 3D Gaussian 
filter in linear rather than logarithmic fc-spacing (note that for low k the smooth power spectrum is 
underestimated). All the lines are divided by the approximation of the no-wiggle power spectrum 
fkpprox given by [34] . 


where P{k) is the initial power spectrum that needs to be smoothed, and Papprox is the initial 
broad band approximation curve (given by e.g. [34, 35]). 

First, we explore a 3d Gaussian filter in linear k spacing: 

where A typically takes a value around 0.05 Mpc//i. For the no-wiggle power spectra we then 
have 

Pnw(fc) = J d^q P(g)PG(|k - q|) 

= + '=')) (A.3) 

The problem that emerges when implementing this method is at the low k limit of the spectra. 
Since the smoothing radius is about the size of wiggle wavelength and since the first wiggle 
starts at the distance closer to the k = 0 than the typical wiggle size, smoothing introduces a 
mismatch of power amplitude at the k = 0 limit (see figure 9). Note also that allowing A to 
vary with scale, A —>• X{k), so that A —?■ 0 as A: —)• 0, is also not an optimal solution since this 
procedure could itself introduce some spurious wiggle-like features, and thus requires more 
complicated filter ansatz. 

Transforming the power spectra to the logarithmic variable, so that the range covers the 
whole real axes, enables us to use simple ID filters. In this case, we get a simple expression 

Pnw J dqiog P (10®°s) exp , (A.4) 
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Figure 10. Grid of B-spline smoothing curves (colored solid lines), varying in spline degree p and the 
number of knots n. Linear power spectrum is shown in dotted back line. For comparison smoothing 
result obtained from applying the ID Gaussian filter is shown as a solid blue line. All the lines are 
divided by the approximation no-wiggle power spectrum Papprox given by [34]. 


where k\og and giog are variables that take values from the whole real axis, and is A a parameter 
with typical value 0.25/i/Mpc. Note that since wiggles are now not equidistantly spaced 
(being more compressed at higher /c), it is beneficial to introduce a slight scale dependence 
of A, increasing the value towards higher k. Result of applying these filters on the linear 
power spectrum is also shown in figure 9. We see that at the low k smooth power spectrum 
is underestimated if one uses the 3D filters in linear spacing. 

Alternative to using the explicit filters to smooth the BAO wiggles is using the basis 
splines (B-splines)^. B-splines offer a convenient way to approximate our wiggle curves. A B- 
spline is a generalisation of the Bezier curve. We can define a knot vector T = [to, ti,, tm], 
where T is a nondecreasing sequence with tj in [0,1], and the set of control points Cq, ..., Cn- 
Degree of the B-spline is then given asp = m — n — 1 and the knots tp+i, ..., tm-p-i are 
called internal knots. The basis functions are defined as 






I 1 if ti < t < ti+i 

I 0 otherwise 
t-ti 


ti-t-j tj 




(A.5) 


where j = 1, 2, ..., p. We can define the B-spline power spectrum approximation as 


Pn,pik) = 

i=0 


(A.6) 


^see e.g. http://mathworld.wolfram.com/B-Spline.html 
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Figure 11. In upper panels the Zel’dovich power spectrum is shown as the blue line, and the resulting 
smoothed version as the orange line. Smoothed version is obtained by computing the Zel’dovich power 
spectrum using the smoothed linear spectrum form Eq. (A.8). Dotted line is linear power spectrum, 
and all the lines are divided by the smooth version of linear spectrum given by Eq. (A.8). In lower 
panels the Zel’dovich power spectrum is shown divided by the smooth version (solid orange line). For 
the reference, the linear power spectrum is also shown divided by the smooth version of itself (black 
dotted line). Left and right panels differ in fc—spacing: linear vs. logarithmic, in order to stress the 
agreement on both large and small scales. 


Knot values determine the extent of the control of control points. In order to achieve smooth¬ 
ing Ci can be determined by the some regression methods and we use a simple linear model 
fit. 

One can construct a family of curves that approximate the original wiggle spectrum. The 
family of such curves is shown in figure 10 for spline degrees p = 2 .. .5 with a various number 
of knots. This family allows us now to impose certain constraints on the final no-wiggle 
spectra Pnw One such constraint can be velocity dispersion 


2 I f d^q Pik) 
'^~3j (27r)3 


(A.7) 


A similar constraint is the erg defined in Eq. (3.3). Note that we can impose several conditions 
at the same time as long the constructed B-spline curves are distributed widely enough. We 
can thus construct the final smooth spectra as a weighted average over the spline curves 


Pnwik) 


^^WsPsik), 


(A.8) 


where summation in s includes all the eligible approximations given by the knot vector T and 
degree of the B-splines p. Imposing the integral constraints like and it| on the Pnw places 
constraints on the weights Wi. 

As a simple example of the smoothing procedure, we can have a look a the Zel’dovich 
power spectrum. In figure 11 we show the results of the Zel’dovich power spectrum computed 
using the power spectrum with and without the wiggles. Note that in this case it is very 
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useful that the smoothed linear spectrum satishes the integral constraints and have the same 

as and a^ values. These ensures that on the large as well as small scales away from the BAO 

region two spectra will agree precisely. 
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